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The spin-orbit coupling systems with a zero magnetic field is studied under the equilibrium sit- 
uation, i.e., without a voltage bias. A persistent spin current is predicted to exist under most 
circumstances, although the persistent charge current and the spin accumulation are identically 
zero. In particular, a two-dimensional quantum wire is investigated in detail. Surprisingly, a per- 
sistent spin current is found to flow along the confined direction, due to the spin precession in 
accompany with the particle motion. This provides an interesting example of constant spin flowing 
without inducing a spin accumulation, contrary to common intuition. 
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In last two decades, the persistent charge currents in 
mesoscopic rings threaded by a magnetic flux have been 
extensively investigated. |ij|2|. L| The persistent current is 
a pure quantum effect appeared in the equilibrium situ- 
ation. In an ideal case it can last forever without dissi- 
pation. In recent years, another subject, the spin-orbit 
(SO) interaction in semiconductors, has also attracted a 
great deal of attention. The SO coupling is an intrin- 
sic interaction, originated from the relativistic effect. By 
using the SO coupling, the electron spins can be conve- 
niently manipulated and controlled experimentally, this 
may have large potential applications in new generation 
of spin-based electronic devices. Due to the SO cou- 
pling, a moving spin is under an equivalent magnetic field 
with the field strength depending on the momentum of 
the particle. Thus, it is nature to ask whether there exists 
a persistent current in a SO coupling system, similar as 
for a system with an external magnetic field. Some recent 
works have begun to explore in this direction. 0, IE 

In this Letter, our goal is to systematically address 
the question: Is there a persistent spin or charge cur- 
rent in a SO coupling system without an external elec- 
tric or a magnetic field? Furthermore, we focus on the 
physical properties of the persistent current if they do 
exist. Before studying a concrete system, we first pro- 
vide a general discussion from symmetry point of view. 
First, consider a system without a magnetic field nor the 
SO coupling. Thus, the Hamilton is H = 2— + V(r), 
with p being the momentum operator and V(r) the po- 
tential energy. This Hamilton contains two symmetries: 
(i) The time-reversal invariance, i.e. [T, H] — with T 
being the time-reversal operator. This symmetry leads 
to the Kramer degeneracy: i.e. if ^(r) is an eigenstate 
of H, the state T^>(r) is also an eigenstate with the 
same energy and < ^\T^f >= 0. 8] (ii) The Hamilton 
commutes with the three Pauli matrices cr, (i = x,y,z), 
i.e. [ui,H] = 0. This symmetry leads to the degener- 
acy: 'J = (</?(r),0) T and = (0, <p(r)) T have the same 
eigenenergy and < >= 0. In the equilibrium case 



with the degenerated energy eigenstates having the same 
occupational probability, it is known that the persistent 
charge current is identically zero because of the time- 
reversal symmetry (i) and so is the persistent spin cur- 
rent because of the symmetry (ii). Second, consider a 
mesoscopic ring threaded by a magnetic flux In this 
system, the time-reversal symmetry (i) is broken but the 
symmetry (ii) survives. Then the persistent spin current 
must vanish because of the symmetry (ii), but the per- 
sistent charge current is in general nonzero. In fact, this 
kind of systems has been extensively investigated and is 
well understood by now. 

011 Finally, consider a sys- 
tem having SO coupling but without a magnetic field. 
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The Hamilton is H — 2— + V(r) + H so . In the present 
system, the symmetry (ii) is destroyed, but the time- 
reversal symmetry still holds since [T,H so ] = 0.0 Hence 
the Kramer degeneracy of VP and T\P exists. These two 
states have exactly opposite charge current densities but 
identical spin current densities. So the persistent charge 
current and the local spin accumulation must be zero ev- 
erywhere in equilibrium. However, the persistent spin 
current usually exists because of the absence of the sym- 
metry (ii). 

Next we consider a concrete system having a SO cou- 
pling. We focus on whether the spin current is indeed 
nonzero in equilibrium, as well as study its physical prop- 
erty if it does exist. Our system is a two-dimensional 
electron gas (2DEG) with the Rashba SO interaction, 
and its Hamilton is: 

H=^^+V(x,y) + ^.(aXp). (I) 
2m n 

The Rashba SO interaction originates from an asym- 
metrical interface electric field in the z-directionjp, llOj 
and its strength a can be varied experimentally. [Tl| The 
potential energy V(x,y) = while < y < L and 
V(x, y) — oo otherwise, i.e. we consider a quasi one- 
dimensional quantum wire as shown in Fig. la. The 
present Hamiltonian can be exactly solved. Due to the 
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fact that [p x , H] = 0, k x is still a good quantum number. 
For a fixed k x , the general solutions of the differential 
equation H^{x : y) = E^>(x,y) are: 



»={1, 2,3,4} 



(2) 



where are constants to be determined by boundary 
conditions, and ^>i(x,y) are: 



^ 4 ( Xj2/ ) = ll^ e -ik-y { ~(K -ik x )/k. 
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kft zkkji, and kR = 



am/h . Using the boundary conditions 0/L) = 0, 
the eigenenergy E can be solved. Fig. 2 shows E versus 
k x and it exhibits a series of sub-band structures due to 
the confinement in the transverse direction. For a nar- 
row quantum wire (e.g. L = 50nm, see Fig. 2a), the sub- 
bands are well apart and separated from each other. The 
curves of E versus k x are of parabolic shape. On the other 
hand, for a wide quantum wire with L = 200nm (see 
Fig. 2b), the sub-bands are very dense, and there are obvi- 
ous mixtures between different sub-bands. These results 
are similar as in the prev ious studies with a parabolic 
confinement potential. [12| 

For given eigenenergy E and k x , i.e., corresponding 
to each point on the curves in Fig. 2, the eigen wave- 
function ^>(x, y) (i.e. the coefficient ai in Eq.(2)) can be 
obtained using the boundary conditions and the normal- 
ization equation. For convenience, we take the state den- 
sity normalization: < >= l/(dE/dk) oc p(E). Then 
the (linear) spin current density j s (x, y) (j s = Re^vs^), 
the angular spin current density (or spin torque) j u (x, y) 
(J u = Re^(ds/dt)^f), and the local spin density s(x,y) 
(s = ^sM^) for every state ^(x,y) can also be calcu- 
lated. Here j s and j u describe respectively the trans- 
lational and the rotational motion of a spin, and they 
satisfythe spin continuity equation jfes = —V • j s + j u 
(see [l3| for a detail discussion). As the last step, sum- 
ming over all the occupied states, the persistent spin cur- 
rent density and the local spin density in the equilib- 
rium can be obtained straightforwardly: e.g. jj(x,y) — 
J dE J2 n f(E)isn(x, V, E) where f(E) is the Fermi distri- 
bution function and n is the sub-band index. 

Next we numerically investigate the linear and angu- 
lar persistent spin current densities and jj. It clearly 
shows that jj and j£ indeed exist in (Rashba) SO cou- 
pling system in the absence of a magnetic field (See 



Fig. 3). We mention again that in the present system 
the persistent charge current and the spin accumulation 
are identically zero because of the time-reversal symme- 
try. Our numerical calculations also confirmed this re- 
sult. The non-zero elements of and j£ are jj yx , jf xy , 
jfxzt and jux- Those non-zero elements are all indepen- 
dent of x because of the translational invariance along 
the x direction. Fig. 3 shows jf yx> jf tXy , and jf xz versus 
y, in which jj iyx and j^ X y are symmetric and jj xz are 
anti-symmetric. j^ x is not shown here because it is re- 
lated to jf yx through the spin continuity equation for a 
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Surprisingly, j* turns out 



to be non-zero. This means that there exists a persistent 
spin current flowing along the y-direction (confined di- 
rection) with the spin pointing in the x-direction. This 
is unexpected since common belief is that if there is spin 
current in the confined direction, there will be a spin ac- 
cumulation at sample boundaries. However, in this sys- 
tem there is no spin accumulation anywhere, but a spin 
current continuously and persistently flows along the y- 
direction. In the following, we provide an understanding 
why this is possible. 

The upshot of the story is that because of the role of 
the angular spin current j£ in the spin continuity equa- 
tion, the spin will precess in accompany with the particle 
motion. As a consequence, the spin current can flow in 
the confined direction without causing spin accumula- 
tion. 

In order to understand the origin of the non-zero jj yx 
in more detail, we investigate the spin current density 
of a single state, namely, the state A in Fig. 2a. Here 
the non-zero elements of j s and j w are still j s<yx , j s , X y> 
js,xz, and j ux . The others elements as well as s x are zero 
for all (x,y). Fig.4a,b show j Siyx , jux, s y and s z of the 
state A. The curves for j s ,xy and j s<xz are similar to s y 
and s z , respectively, so they are not shown here. It is 
surprising that there is a non-zero spin current flowing 
along the y-axis with spin pointing to the a;-direction, 
but the element s x is exactly zero everywhere. 

Before to explain the spin current density, let us first 
study what is the electron motion in the state A. With- 
out the Rashba coupling, i.e., kn — > 0, its wave-function 
is V(x,y) = Ne ik * x sin{ny/L)(i,l) T = K^x^wy/L _ 
e-™ y / L )(i, 1) T = * + + This is a state with elec- 
tron spin towards — y direction, propagating along the x 
axis while reflecting back and forth in the transverse y- 
direction, as shown in Fig. la. Here \&+ and represent 
the traveling waves along +y and — y directions, respec- 
tively. While fee ^ 0, the wave-function "J can still be 
separated into and with fy + = +03^/3 and 
\I r _ = 02^2 + 0,4^4- So the electron motion looks still 
like the case of kn = 0. However, its spin is no longer 
fixed at — y direction anymore, and the electron spin pre- 
cesses while an electron makes its motion because of the 
existence of the Rashba SO coupling, Let us describe 
the spin precession in detail for the state A. Assume an 
electron at the beginning is on the edge of y = mov- 
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ing in the +y direction (see the point 1 in Fig. la). This 
state is described by wave-function \I/+(x, 0). The spin 
of \f r +(x, 0) is mainly pointing to —y with a slight tilt 
towards +z (see the arrow 1 in Fig. lb). When an elec- 
tron makes the motion from 1 — * 2 — * 3 (described by 
^f+(x,y)), then 3 — > 4 — > 5 (described by \t_(x, y)), the 
spin vector makes the rotation as shown in Fig. lb. When 
the electron comes back to the original side, the spin also 
precesses back to the original direction it started with. 
The above process repeats itself. 

With the above knowledge, we are ready to explain the 
origin of the spin current. We still use the state A as an 
example. Note that everywhere inside the sample one can 
split the wave- function: \t(x, y) = $f+(x,y) + \I>_(x, y). 
For the +y-direction moving wave its spin element 
s x is negative (e.g. see the arrow 2 in Fig. lb), so it in- 
duces a negative j s ,yx- For the — y-direction moving wave 
its spin element s x is positive (e.g. see the arrow 4 in 
Fig. lb), so it also has a negative j s , y x because the motion 
is along —y. Therefore, a negative j s ,yx emerges for the 
state A (see Fig. 4a). The behavior of other elements of 
js(x, y), ju>{x, y), and s(x, y) can also be explained in the 
same manner. Here we emphasize two points: (i) The be- 
haviors of all elements of j s , j u , and s can be understood 
with the above picture of the electron motion accompa- 
nied by its spin precession, (ii) Due to the fact that the 
motion is periodic and repeated, this motion does not 
change the local spin density, and the spin current can 
persistently exist. 

In the above analysis, the electron is in the lowest sub- 
band. The electron motion in a higher sub-band is simi- 
lar, only its spin precession is in a more complicated man- 
ner. The electron motion plus the spin procession can 
also generate nonzero j StVX , j StXy , j StXZ , and j ux . Fig.4 
(c) and (d) show j S)VX and j ux for the state E (see Fig. 2a) 
in the 2nd sub-band and the state F (see Fig. 2b) in the 
4th sub-band, respectively. We emphasize that those re- 
sults and the nature of the electron motion show that the 
wave-functions are extended in the transverse direction 
and they do not form edge states. This is quite different 
from the case with a strong magnetic field, in which edge 
states appear. 

The state B (see Fig. 2a), the other state in the 1st 
sub-band, describes a similar electron motion (shown in 
Fig. la, c), but the spin is mainly towards +y axis and 
it tilts to — z axis at beginning (point 1). So that its 
spin current elements are opposite in sign with those of 
the state A, including j s<yx and j^x which are shown in 
Fig. 5a, b. Thus, they partially cancel each other out, but 
the remaining net spin currents are still quite large. Also 



note that there are total four degenerated states (A, B, 
C, and D) at E = 0.005 and L = 50nm. The other two 
states C and D are the time- reversal states of the states A 
and B, respectively. Therefore, the spin currents of C and 
D are exactly the same with those of A and B, but the 
local spin densities and charge currents are exact opposite 
with those of A and B. So that the total spin current of 
those four states doubles that of A and B, and the local 
spin density or the charge current vanishes. While at a 
higher E or for a wider sample width L, more sub-bands 
are involved and more degenerated states emerge. Fig. 5c 
shows the spin currents j s ,yx and j ux for E = 0.005 and 
L = 200nm with 4 active sub-bands. The total spin 
current of those states still remain a large value. 

Integrating over energy E, the persistent spin currents 
are obtained, which are shown in Fig. 3. From the above 
analysis, the spin current j s ,yx flowing in the confined di- 
rection can indeed exist constantly and it does not induce 
any spin accumulation. This is due to the spin precession 
during the electron motion, which makes a nonzero angu- 
lar spin current j ux . Finally, we make several remarks, 
(i) All electrons in the system, not only the ones near 
the Fermi surface, contribute to the persistent spin cur- 
rent. This is quite different from a voltage-biased system, 
in which only the electrons near the Fermi surface con- 
tribute to a current, (ii) It is a universal feature that a 
persistent spin current appears in the (Rashba) SO cou- 
pling systems. For example, in a mesoscopic ring or a 
disc device, j s and j u are usually nonzero. Moreover, 
impurities in the device will not destroy the persistent 
spin current, the spin flow will only be modified by the 
impurity potentials, (iii) The persistent spin current can 
induce an electric field, which offer a way of detecting the 
persistent spin current. [l3lll5| 

In summary, we predict that a persistent spin current 
will commonly appear in a spin-orbit coupling sample in 
the absence of a magnetic field under zero voltage bias, 
although the spin accumulation and the persistent charge 
current are identically zero. Moreover, a spin current j s is 
found to flow along the confined direction in the Rashba 
SO coupling two-dimension quantum wires. This implies 
that a spin current j s does not necessarily related to spin 
accumulation. 
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FIG. 1: (color online) (a) Schematic diagram for a two- 
dimensional quantum wire device. The arrows describe the 
direction of particle motion for a state with a positive k x . (b) 
and (c) are schematic diagrams for the spin precessions in the 
states A and B in Fig. 2, respectively. 

FIG. 2: (color online) The dispersion relation, E vs. k x , for 
L = 50nm (a) and 200nm (b) with kit = 0.01/nm. 

FIG. 3: (color online) jJ tXy (blue dashed line), jf xx (red dot- 
ted line), and jj yx (black solid line) vs y for L — 50nm and 
E F = 0.02 (a), 'and £ = 200nm and E F = 0.005 (b), re- 
spectively. The other parameters are ku = 0.01/nm and the 
temperature T = 0. Here p = 2mEF/h 2 is the electron area 
density. 

FIG. 4: (a) and (b) show j s , V x (solid line in (a)), j wx (dotted 
line in (a)), s y (solid line in (b)), and s z (dotted line in (b)) 
vs. y for the state A in Fig. 2a. (c) and (d) show j s , yx (solid 
line) and j ux (dotted line) vs. y for the state E in Fig. 2a and 
the state F in Fig. 2b, respectively. 

FIG. 5: (color online) Show j s , yx (a) and j wx (b) vs. y for the 
state A (black solid line), the state B (black dotted line), and 
the sum of the two (red solid line), (c) Show the total j 3 , yx 
(black solid line) and j ux (red dotted line) of the 16 states for 
E = 0.005 and L = 200nm (see Fig.2b) vs. y. 
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